Magnetic order in GdBiPt studied by x-ray resonant magnetic scattering 
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Rare earth (K) half-Heusler compounds, iiBiPt, exhibit a wide spectrum of novel ground states. Recently, 
GdBiPt has been proposed as a potential antiferromagnetic topological insulator (AFTI). We have employed x- 
ray resonant magnetic scattering to elucidate the microscopic details of the magnetic structure in GdBiPt below 
Tjv = 8.5 K. Experiments at the Gd L2 absorption edge show that the Gd moments order in an antiferromagnetic 
stacking along the cubic diagonal [111] direction satisfying the requirement for an AFTI, where both time- 
reversal symmetry and lattice translational symmetry are broken, but their product is conserved. 

PACS numbers: 75.25.-j, 75.50.Ee, 73.20.-r 



The discovery of three-dimensional topological-insulating 
states in binary alloys (Bii_a;Sba;)i^ and compounds (Bi2Se3, 
Bi2Te3, Sb2Te3),-'"^ which feature an insulating gap in the 
bulk but with topologically protected conducting states on 
the surfaces or edges, has opened a new frontier for funda- 
mental condensed matter physics research.^ As pointed out 
in several papers, the novel properties of this class of mate- 
rials offer potential for technological breakthroughs in quan- 
tum computing and magneto-electronic applications.'^'* Over 
the past year, attention has turned towards investigations of 
new phenomena that arise when topological insulators (TI) 
also manifest, or are in close proximity to, other phenom- 
ena including magnetic order and superconductivity.^^— Re- 
cently, the Heusler and half-Heusler compounds have been 
subject to intense scrutiny because of their potential as TI with 
tunable electronic properties.""''* Specifically, Mong et at.— 
have proposed that GdBiPt may provide the first realization 
of an antiferromagnetic topological insulator (AFTI), where 
both time-reversal symmetry and lattice translational symme- 
try are broken, but their product is conserved. Predictions 
for this class of TI include gapped states on some surfaces, 
gapless states on others, and novel one-dimensional metallic 
states along step edges on the gapped surfaces.'" 

More generally, Heusler and half-Heusler compounds ex- 
hibit a wide spectrum of novel ground states.'^ The rare 
earth (R) half-Heusler compounds, i?BiPt, feature magnetic 
ordering (GdBiPt),'^ superconductivity (LaBiPt, YBiPt)'^ '^ 
and heavy-fermion behavior (YbBiPt)'^ . Although the low- 
temperature ground states of the i?BiPt system (for R = Ce, 
Nd, Sm, Gd, Tb, Dy, Ho, Er, Tm, and Yb) have been char- 
acterized as antiferromagnetic through thermodynamic and 
transport measurements, there have been few magnetic struc- 
ture determinations for this series.^" GdBiPt has the highest 
T/v of the series at approximately 8.5 K'^ and, since the or- 
bital angular moment L = for the 5-state Gd ion, the mag- 
netic structure in the absence of crystalline electric field ef- 
fects may be directly investigated. However, the high neutron- 
absorption cross section for naturally occurring Gd is prob- 
lematic for conventional magnetic diffraction experiments. 

Here we describe the magnetic order of GdBiPt below 
Tjv = 8.5 K determined by x-ray resonant magnetic scattering 
(XRMS) at the Gd L2 absorption edge. GdBiPt crystallizes 




FIG. 1 : (Color online) Crystal structure of GdBiPt. 



in the MgAgAs-type structure (cubic space group F 4 3 m, 
a = 6.68 A with Gd, Bi and Pt at the 4c, 4d and 4a sites, re- 
spectively; see Fig. [T]ii2ii2^ The structure may be viewed as 
three sets of elementally pure, interpenetrating face-centered 
cubic lattices. We find that the commensurate magnetic or- 
der doubles the cubic unit cell along the diagonal [1 1 1] di- 
rection, characterized by a propagation vector q„j = (5^5)^ 
so that alternating ferromagnetic (111) planes of Gd are an- 
tiferromagnetically coupled along the [111] direction. This 
structure is quite similar to the model B magnetic structure 
for an AFTI via spin-orbit coupling as described by Mong et 
fl/.,'** but we find that the moment direction in GdBiPt is not 
parallel to the magnetic propagation vector as is found, for 
example, in MnSbCu^^ or CeBiPt.^o 

Single crystals of GdBiPt were solution-grown using a 
Bi flux and emerged with sizeable facets perpendicular to 
the [001] direction and smaller facets perpendicular to 
[1 1 1]. High-purity Gd (obtained from Ames Laboratory), 
Pt, and Bi were placed in an alumina crucible in the ratio 
Gd : Pt : Bi = 3 : 3 : 94, sealed in a silica ampule, and slowly 
cooled from 1 170°C to 600°C over 200 hours. At 600°C, the 
excess Bi solution was decanted from the GdBiPt crystals.— 
The dimensions of the single crystal studied in the XRMS 
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FIG. 2: (Color online) Resonant magnetic scattering from the GdBiPt 
single crystal, (a) Rocking scans (9) through the (i -i ^) magnetic 
peak position above (open circles) and below (filled circles) Tat taken 
in CT-TT scattering geometry, (b) Energy scan through the Gd L2 ab- 
sorption edge at the (i -i ^) magnetic peak position at T = 4.7 K 
(blue filled circles) along with the measured x-ray fluorescence from 
the sample (black filled circles). 
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FIG. 3: (Color online) (a) M/ H and its temperature derivative for 
GdBiPt. (b) The magnetic intensity measured while scanning tem- 
perature at the maximum of the (5 ^ |) diffraction peak without a 
polarization analyzer (open small circles) and the integrated intensity 
of the (^ -| ^) diffraction peak measured at selected temperatures 
using the polarization analyzer (open large squares). The solid line 
is a power law fit to the integrated intensity data as described in the 
text. 



measurements were approximately 3x3x2 mm^ with a large 
as-grown facet perpendicular to [GDI]. The measured mo- 
saicity of the crystal was less than 0.01 degrees full-width-at- 
half-maximum (FWHM), attesting to the high quality of the 
sample. The XRMS experiment was performed on the 6ID-B 
beamline at the Advanced Photon Source at the Gd L2-edge 
{E = 7.934keV). The incident radiation was linearly polarized 
perpendicular to the vertical scattering plane ((T-polarized) 
with a beam size of 0.5 mm (horizontal) x 0.2 mm (vertical). 
In this configuration, dipole resonant magnetic scattering ro- 
tates the plane of linear polarization into the scattering plane 
(TT-polarization). For some of the measurements, pyrolytic 
graphite PG (0 6) was used as a polarization analyzer to sup- 
press the charge and fluorescence background relative to the 
magnetic scattering signal. For measurements of the magnetic 
reflections, the sample was mounted at the end of the cold fin- 
ger of a closed-cycle cryogenic refrigerator with the {HHL) 
plane coincident with the scattering plane. 

Measurements of the diffraction from the sample performed 
in the cr-vr scattering geometry using the PG (0 6) polariza- 
tion analyzer are shown in Fig. |2] For temperatures above 
T/v = 8.5 K, only Bragg peaks consistent with the chemical 
structure^ii^ of GdBiPt were observed. However, upon cool- 



ing below Tjsi, additional Bragg scattering at half-integer val- 
ues of (HKL) was found as shown in Fig.|2|a). The magnetic 
origin of these peaks was confirmed by energy scans through 
the Gd L2 absorption edge and from the temperature depen- 
dence of the diffraction peak intensity as described below. 

The energy scan Fig.|2|b) was performed with the diffrac- 
tometer set at the magnetic peak position and is typical of 
resonant magnetic scattering at the L edges of rare-earth 
compounds. At the L2 edge of rare-earth elements, the reso- 
nance primarily involves electric dipole {E\) transitions from 
the 2pi core level to the empty 5d states, seen as the strong 
line just at, or slightly below the maximum in the measured 
fluorescence intensity. The weaker feature below the E\ res- 
onance in Fig. |2b) is likely due to the electric quadrupole 
(E2) transition from the 2pi core level to the 4/ states that 
are pulled below the Fermi energy because of the presence of 
the core hole in the resonance process. 

The temperature dependence of the magnetic scattering, 
along with the corresponding magnetization measurements 
performed on a sample from the same batch using a Quantum 
Design Magnetic Properties Measurement System, are shown 
in Fig. [3] The magnetic order parameter was measured at the 
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FIG. 4: (Color online) Integrated intensity in azimuth scans through 
the (-i -i ^) magnetic Bragg peak. Measured data are depicted 
by full black circles. Full and dashed lines represent calculations for 
selected magnetic moment directions in a single magnetic domain, 
and for intensity from equally populated domains averaged over the 
three possible symmetry-equivalent magnetic moment orientations, 
respectively. 



(i i |) peak position as the sample temperature was increased 
during a temperature scan in the absence of the polarization 
analyzer. These data were supplemented by measurements of 
the integrated intensity of the (i -i ^) magnetic Bragg peak 
at selected temperatures and with polarization analysis. The 
line in Fig.[3tb) describes a fit to the integrated intensity data 
using a power law of the form / ^ (1 — yielding 
T/v = 8.52±0.05K and /3 = 0.33±0.02. The close proximity 
of Tat determined from our scattering measurements and the 
peak in d[(M/i/)T]/dr (Ref.26) at T = 8.6 K, again con- 
firms the magnetic origin of the Bragg scattering with a prop- 
agation vector of = (i I i). Systematic M versus H mea- 
surements (not shown) demonstrate in addition that no spon- 
taneous ferromagnetic moment is present. 

Having established the nature of the magnetic ordering in 
GdBiPt, we now describe our attempt to determine the direc- 
tion of the ordered magnetic moment. The angular depen- 
dence of the resonant magnetic intensity /(i/') for the incident 
CT-polarized beam depends upon the component of the mag- 
netic moment along the scattered beam direction and can be 

written as /(V')(q, a, /3) = C [m ■ k'(V')(Q)]^ ^(V')(q, a, /3) 
where C is an overall scale factor that accounts for the res- 
onant scattering matrix element and incident beam intensity, 
m and k' represent the magnetic moment and scattered beam 
directions, respectively, and A accounts for the absorption 
correction.-^ The sample geometry required off-specular scat- 



tering measurements of the magnetic peaks. That is, the angle 
a of the incident beam k with respect to the sample surface 
is different from the angle /3 of the outgoing beam k' with re- 
spect to the sample surface.^** For the azimuth angle i\} scans 
shown in Fig. H] the diffractometer was set at the position of 
the magnetic Bragg peak and the crystal was rotated around 

the scattering vector Q = k' — k thereby rotating k' with 
respect to m while leaving Q fixed. This yields an azimuth 
dependence of the intensity which is specific to a given mag- 
netic moment direction. Note, that the absorption correction 
A also depends on the azimuth angle V'. 

For a cubic lattice, the determination of the ordered mo- 
ment direction is significantly complicated by the presence of 
domains that arise from symmetry-equivalent magnetic propa- 
gation vectors and moment directions. For the observed mag- 
netic Bragg peaks at (iJ K L) with H, K, and L half integers, 
four symmetry-equivalent | 5} propagation vectors exist: 
(i i i), (-i -i i), (i -i -i), and (-!_ \ -i). Fortunately, 
for GdBiPt in the cubic space group F 4 3 m, only one prop- 
agation vector contributes to a particular magnetic reflection 
[e. g. the magnetic Bragg peak (-i -i -^) is generated by the 
propagation vector q„j = (-\ -\ i) from the (006) zone cen- 
ter]. The measured data for the {-\-\ ^) magnetic Bragg 
peak show two important features in the azimuth scan pre- 
sented in Fig.m a distinct minimum with almost no intensity 
close to = 0, and an increase in intensity by more than an 
order of magnitude as t/j is varied by ±30°. Both features 
are in strong contrast to the expected dependence of the in- 
tensity for magnetic moments parallel to the propagation vec- 
tor = (-i -i i) as illustrated in Fig. |4]by the bold black 
line with a maximum close to = 0. Therefore, we can ex- 
clude that the moments are parallel to the propagation vector 
in GdBiPt. 

In Fig. m calculated curves are also shown for other mo- 
ment directions. Unfortunately, for each of the depicted mo- 
ment directions, three different symmetry-equivalent orienta- 
tions can occur yielding three magnetic domains. The dashed 
lines in Fig. |4] represent the calculated dependence of the 
intensity if we include all such domains for a given moment 
direction with equal population. We again find poor agree- 
ment between the domain-averaged calculations for moments 
along the set of {1 1 1}, {1 1 0} and {0 1} directions. 
However, calculations assuming the presence of only a sin- 
gle domain within the probed volume, with one specified mo- 
ment direction (either [-1 1 1] or [-1 1 0] for the (-i -i ^) 
Bragg peak in Fig.|4|l come much closer to describing the mea- 
sured data. This behavior clearly indicates that the magnetic 
domains are large; smaller than the footprint of the incident 
beam on the sample (approximately 0.5 x 0.5 mm'^), but of the 
same order of magnitude. Similar large magnetic domains 
have been noted in previous XRMS work on GdNi2Ge2 as 
wellf22 Nevertheless, a unique determination of the moment 
direction is not possible based on the available data. A more 
precise determination of the moment direction may be pos- 
sible from measurements with much smaller incident beam 
dimensions and/or control of domain populations.^' 

Summarising the experimental results, below T/v = 8.5 K 
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the magnetic Gd moments order in a commensurate antifer- 
romagnetic structure in GdBiPt that can be described as dou- 
bhng the cubic unit cell along the diagonal [111] direction, 
so that alternating ferromagnetic (111) planes of Gd are an- 
tiferromagnetically coupled along the [111] direction. The 
moments are not aligned parallel to this diagonal [111] direc- 
tion. 

In contrast to GdBiPt, CeBiPt is an antiferromagnet char- 
acterized by a propagation vector q„j = (100) and the or- 
dered moments are collinear with the propagation vector 
along [1 00],^*^ but with a reduced moment that may, in part, 
be attributed to crystalline electric field (CEF) effects. Un- 
fortunately, XRMS measurements do not allow a direct ex- 
traction of the ordered moment in GdBiPt, but earlier specific 
heat measurements"' estimated an entropy of ^0.8 7?. In 8 as- 
sociated with the magnetic transition close to the value ex- 
pected for full moment ordering without CEF effects. The en- 
tropy associated with the corresponding magnetic transitions 
for the Nd, Tb and Dy compounds were considerably less than 
7?^ln(2 Jh-1) expected for the full Hund's rule J multiplet, in- 
dicating the importance of CEF effects in these compounds. 
The magnetic structures for B = Nd, Sm, Tb, Dy, Ho, Er, Tm, 
and Yb, have not yet been identified by neutron or XRMS 
measurements and such measurments are planned. 

Finally, we comment on our results in light of the proposal 



that GdBiPt may be an AFTI candidate. The AFTI state may 
be derived from either magnetic ordering in a pre-existing 
strong TI (model A in Ref.llOl) or, alternatively, for specific an- 
tiferromagnetic ordering schemes that induce spin-orbit cou- 
pling in the system (model B in Ref. 10). Given previous 
ARPES measurements^ above Tn, which do not find direct 
evidence for band inversion in GdBiPt, it seems unlikely that 
GdBiPt is itself a strong TI. However, the magnetic structure 
determined here is consistent with the alternative model B pre- 
sented by Mong et al. The doubling along the cubic diagonal 
direction represents the broken lattice translational symme- 
try (by order of two) and the ordering of each magnetic mo- 
ment breaks the time-reversal symmetry, however, the product 
of both symmetry operations is conserved for the determined 
magnetic order. 
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